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The non-Hermitian PT-symmetric quantum-mechanical Hamiltonian H = p 2 + 
x 2 (ix) 6 has real, positive, and discrete eigenvalues for all e > 0. These eigenval- 
ues are analytic continuations of the harmonic-oscillator eigenvalues E n = 2n + 1 
^5 ( n = 0> 1, 2, 3, . . .) at e = 0. However, the harmonic oscillator also has nega- 

Q\ tive eigenvalues E n = —2n — 1 (n = 0, 1, 2, 3, . . .), and one may ask whether it 

is equally possible to continue analytically from these eigenvalues. It is shown in 
this paper that for appropriate 'PT-symmetric boundary conditions the Hamilto- 
nian H = p 2 + x 2 (ix) e also has real and negative discrete eigenvalues. The negative 
eigenvalues fall into classes labeled by the integer N (N = 1, 2, 3, . . .). For the iVth 
class of eigenvalues, e lies in the range (4./V — 6)/3 < e < AN — 2. At the low and high 
ends of this range, the eigenvalues are all infinite. At the special intermediate value 
e = 2N — 2 the eigenvalues are the negatives of those of the conventional Hermitian 
Hamiltonian H = p 2 + x . However, when e ^ 2N — 2, there are infinitely many 
complex eigenvalues. Thus, while the positive-spectrum sector of the Hamiltonian 
H = p 2 + x 2 (ix) £ has an unbroken VT symmetry (the eigenvalues are all real), the 
negative-spectrum sector of H = p 2 + x 2 {ix) £ has a broken VT symmetry (only some 
of the eigenvalues are real) . 
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CG I. INTRODUCTION 

In 1993 an observation was made pQ regarding the eigenvalues of the conventional 
quantum-harmonic-oscillator Hamiltonian 

H = p 2 + wV, (1) 
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where a; is a real positive parameter. It was noted that while the standard eigenvalues of H 
are real and positive 

E n = (2n + l)co (n = 0, 1, 2, 3, ...), (2) 

if (fil) and (J2J) are analytically continued in the complex-w plane from positive to negative u>, 
the eigenvalues E n all become negative even though the Hamiltonian (fl| appears to remain 
unchanged. Surprisingly, the harmonic-oscillator Hamiltonian (fl|, which is a sum of squares, 
also possesses an infinite set of negative eigenvalues. 

A careful treatment of the boundary conditions on the eigenfunctions is required to 
explain the appearance of negative eigenvalues. The conventional boundary conditions as- 
sociated with the Schrodinger eigenvalue differential equation 

- ip"(x) + wVfy(z) = Eip(x) (3) 

are that ip(x) — > as x — > ±oo on the real axis. These boundary conditions hold not just on 
the real axis but in Stokes wedges centered about the positive-real and negative-real axes in 
the complex- x plane [2]. Specifically, the eigenfunctions ip(x) vanish exponentially rapidly 
in the Stokes wedges defined by — 7r/4 < argx < 7r/4 and at — 57r/4 < argx < — 37r/4. 
As the parameter u rotates in the complex-cu plane in the positive (anticlockwise) direction 
from positive to negative values, the Stokes wedges in the complex-a; plane rotate by 7r/2 
in the negative (clockwise) direction and end up centered about the positive- and negative- 
imaginary axes. We thus ascertain that the harmonic-oscillator Hamiltonian (fTl) has two sets 
of eigenvalues: positive eigenvalues, which arise when the boundary conditions are imposed 
in a pair of Stokes wedges centered about the real- a: axis, and negative eigenvalues, which 
arise when the boundary conditions are imposed in a pair of Stokes wedges centered about 
the imaginary-x axis. 

The negative eigenvalues can be obtained directly by making the transformation 

x = it. (4) 

This simple transformation replaces E in (J3J) by —E and simultaneously replaces the bound- 
ary conditions on the real- a: axis with boundary conditions on the imaginary-^ axis. 

The question addressed in this paper is whether the negative-eigenvalue problem for the 
harmonic oscillator can be extended into the complex domain in a PT-symmetric fash- 
ion. To understand this question let us recall that in order to construct a "PT-symmetric 
Hamiltonian, one begins with a conventional Hermitian Hamiltonian, such as the quantum- 
harmonic-oscillator Hamiltonian, H = p 2 + x 2 , and then introduces a parameter e to extend 
the Hamiltonian into the complex domain while preserving its VT symmetry. The standard 
example of such a Hamiltonian is [SHE] 

H = p 2 + x 2 (ix) e (ereal). (5) 

As e varies smoothly away from 0, the eigenvalues of this Hamiltonian smoothly deform 
away from the harmonic-oscillator eigenvalues. If we begin with the positive harmonic- 
oscillator eigenvalues, the resulting discrete eigenvalues remain real and positive for all e > 
and the eigenvalues grow with increasing e, as shown in Fig. [T] When the eigenvalues of 
a "PT-symmetric Hamiltonian are real, the Hamiltonian is said to have an unbroken VT 
symmetry. A Hamiltonian with unbroken VT symmetry represents a physically viable and 
realistic quantum system, and such systems have been repeatedly observed and studied in 
laboratory experiments [9H20]. 




FIG. 1: Energy levels of the parametric family of Hamiltonians H = p 2 + x 2 (ix) 6 (ereal). When 
e > 0, the eigenvalues are all real and positive, and increase with increasing s. When e decreases 
below 0, the eigenvalues disappear into the complex plane as complex conjugate pairs. Eventually, 
only one real eigenvalue remains when e is less than about —0.57, and as e approaches —1 from 
above, this eigenvalue becomes infinite. 



What happens if we begin with the negative harmonic-oscillator eigenvalues at e — 0, 
instead of the positive eigenvalues, and smoothly increase or decrease si Do the eigenvalues 
remain real and negative? More generally, what happens if we begin with the negative 
eigenvalues that one obtains when e = 2N — 2 (N = 1, 2, 3, . . .) [see ([6|], and then vary 



e? Do the eigenvalues all remain real and negative? We will see that if we begin with the 
negative-real eigenvalues at e — 2N — 2, the smallest-negative eigenvalue remains real and 
negative but only for a finite range of e and not an infinite range of e: 



AN -6 



< e < AN - 2 



(N 



1, 2, 3, 



(6) 



As e approaches the upper and lower edges of the Nth region, this eigenvalue approaches 
— oo. The larger-negative eigenvalues eventually become complex. When e ^ 2N — 2, there 
are always a finite number of real negative eigenvalues and an infinite number of complex 
eigenvalues. The behaviors of the eigenvalues in the first three regions N — 1, 2, 3 are shown 
in Fig. [2j 

This paper is organized as follows: In Sec. |TT] we explain the role of the Stokes wedges for 
the PT-symmetric negative-eigenvalue problem for the Hamiltonian in (pi). Then in Sec. Ill 
we use WKB to obtain accurate numerical approximations to the real eigenvalues. We show 
that WKB provides a clear explanation for why the eigenvalues in the Nth region diverge 
at the lower and upper ends of the region, namely, that the turning points rotate out of the 
Stokes wedges in which the eigenvalue problem is posed. In Sec. IV we make some brief 
concluding remarks. 
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FIG. 2: Negative eigenvalues of the Hamiltonian H = p 2 + x 2 (ix) 6 (e real) in the first three regions 
of e, —2/3 < e < 2 (dots), 2/3 < e < 6 (squares), and 2 < e < 10 (diamonds), corresponding to 
N = 1, 2, and 3 in (raj). (In the electronic version the N = 1, N = 2, and iV = 3 eigenvalues are 
plotted as blue dots, red squares, and purple diamonds, respectively.) At e = 0, 2, 4 the eigenvalues 
are the exact negatives of the conventional positive eigenvalues of the Hermitian Hamiltonians 
H = p 2 + x 2 , H = p 2 + x 4 , H = p 2 + x 6 . When e 7^ 2iV — 2, there are only a finite number of real 
negative eigenvalues and the remaining eigenvalues are complex. 



II. STOKES WEDGES 

To construct a PT-symmetric extension of the Hamiltonian if in (J5l), we must recall 
the effects of V and T on the complex coordinate x. Under space reflection (parity) the 
coordinate x changes sign, V '■ x — ¥ —x, and under time reversal i changes sign, T : x — > x*. 
Therefore, under the combined VT operation the complex coordinate x is reflected about 
the imaginary axis. Thus, as the Schrodinger eigenvalue problem 



- r(x. 



+ x 2 (ix^ £ 



ip(x) = Eip(x) 



(7) 



is posed in the cut complex- a; plane, the requirement of VT symmetry demands that the 
cut lie on the imaginary-x axis. We take the cut to lie on the positive-imaginary axis (see 
Fig. [3]) because for the usual positive-eigenvalue solutions to the Hamiltonian ^ the cut 
was originally taken to lie on the positive-imaginary axis [3]. 

The eigenvalue problem (171) is then posed on a three-sheeted Riemann surface. As is shown 
in Fig. pi on sheet —1 the complex argument of x lies in the range — \ix < argx < — §7r, on 



sheet (the principal sheet) the range is — |vr < argx < W, and on sheet 1 the range is 



\ix < argx < §7r. 
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FIG. 3: Configuration of Stokes wedges and turning points for the case N = 1. As e increases, 
the right wedge (solid lines), which begins on sheets and 1, rotates clockwise and downward, 
as indicated by the solid double arrow and the wedge lies entirely on sheet when e > 1. The 
wedge becomes thinner as it rotates, and its opening angle vanishes as e — > oo; at this point the 
center of the wedge lies at the angle —n/2. The behavior of the left wedge (dotted lines) mirrors 
the behavior of the right wedge, and its rotation is indicated by the double dotted arrow. The 
eigenvalue problem associated with these wedges has real eigenvalues when the turning points lie 
inside the wedges, and this occurs only when —2/3 < e < 2. 

Using the techniques described in detail in Ref. [3], we identify the Stokes wedges in 
which the boundary conditions of the differential-equation eigenvalue problem are imposed. 
The centers of the right and left wedges are 
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(8) 



These angles lie on the principal sheet and are "PT-symmetric reflections of one another. 
The opening angle of the wedges is given by 
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The angular locations of the upper and lower edges of the right wedge are 
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and the angular locations of the upper and lower edges of the left wedge are 
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(9) 
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Note that the right wedge extends onto sheet 1 and the left wedge extends onto sheet —1 until 
e becomes larger than 1. The wedges become thinner and rotate downward as e increases, 
and as e — > oo, the two wedges become infinitely thin and approach — §7r. 

For positive eigenvalues E, the turning points satisfy the equation x 2 (ix) e = E and the 
turning points lie on the real axis when e = 0. However, for negative eigenvalues E = —\E\, 
the turning points satisfy the equation 

{txf+z = -\E\. (12) 

There is a pair of turning points located at the angles 

2 1 2 1 

"turning point— R ~Z . <f o ' "turning point— L ~Z . <! o V / 

Thus, when e — 0, the turning points lie on the imaginary axis and sit in the center of each 
wedge. 

As £ varies, both the wedges and the turning points rotate in the complex- x plane, but 



the turning points rotate faster than the wedges. Thus, the turning points in (13) only lie 
inside of the wedges for the range of e 

-\<e<2. (14) 



Therefore, as we explain in Sec. Ill, the WKB asymptotic estimate of the eigenvalues breaks 
down when e is not in this range. This phenomenon of turning points entering and leaving 
wedges is discussed in Ref. [21]; this phenomenon does not occur for the case of the positive 
eigenvalues discussed in Ref. [3]. 



There are infinitely many solutions to the turning-point equation (12). The angular 



distance between successive turning points is ^j-vr. The integer N labels the turning points 
and the iVth pair of turning points lies at the angles 

^turning points = ± „ ^ ~ 9^ (N — 1, 2, 3, . . .). (15) 

The iVth turning point lies in the Nth pair of wedges given by 

2j\T + 3 1 

Supper edges = ^^fT ^ ~~ O^ (N = 1, 2, 3, . . .), 

Slower edges = ± ^~ n ~ ^ ( N = 1 > 2 ' 3 ' • • ')• ( 16 ) 

The condition that the iVth turning point lies in the iVth pair of wedges is (J6J). [The region 



of £ in (14) corresponds to N — 1. All the turning points and wedges collapse to the angle 



— ^tt as £ — > oo. 



III. WKB CALCULATION OF THE NEGATIVE EIGENVALUES 

We can use WKB to obtain an approximate formula for the nth eigenvalue in the iVtli 
range of e. To do so, we simply evaluate the WKB quantization formula 

/ dx y/E n - V(x) ~ in + 1/2)tt (n ->■ oo), (17) 

J X1 
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FIG. 4: Comparison of the negative eigenvalues in the N = 1 region of e, —2/3 < e < 2, with the 



WKB asymptotic formula ( 18 ) for the eigenvalues. The eigenvalues are plotted as dots (colored blue 



in the electronic version) and the WKB formula is plotted as solid curves (colored light blue in the 
electronic version). At e = the eigenvalues are exactly —1, —3, —5, ..., which are the negatives 
of the conventional positive eigenvalues of the harmonic-oscillator Hamiltonian H = p 2 + x 2 . Note 
that when e^O there are only a finite number of negative eigenvalues. 



where the path of integration runs from the left turning point to the right turning point 
The result is 

-(«+i)r(i + ^r :!+ -'-' 1 " 



E n 



r (i + M - (ff ) 



(18) 



Note that this formula breaks down when e approaches the lower and upper endpoints in 
(|6j) because the cosine vanishes at these points. 

We have computed the eigenvalues in the Nth region (N = 1, 2, 3) numerically by 
integrating along the centers of the wedges and matching at the origin. This matching 
requires that the path go around the branch point at the origin, and thus the procedure is 
reminiscent of the toboggan contours studied by Znojil [22]. The numerical values of the 
eigenvalues are compared with the WKB prediction in (18) in Fig. [4] for N — 1, Fig. [5] for 
N = 2, and Fig. [6] for N = 3. Note that WKB is most accurate when the turning points 
lie exactly in the centers of the wedges. When the turning points do not lie in the centers 
of the wedges, the accuracy of the WKB approximation at first increases with increasing n, 
but the accuracy eventually decreases and the WKB approximation fails entirely when the 
eigenvalues become complex. 
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FIG. 5: Comparison of the negative eigenvalues in the N = 2 region of e, 2/3 < e < 6, with the 
WKB asymptotic formula for these eigenvalues. The eigenvalues are plotted as squares (colored red 
in the electronic version) and the WKB formula is plotted as solid curves (green in the electronic 
version) . At e = 2 the eigenvalues are exactly the negatives of the conventional quartic anharmonic- 
oscillator Hamiltonian H = p 2 + x 4 . When e / 2 there are only a finite number of negative 
eigenvalues. The snake-like behavior of the eigenvalues as functions of e is similar to what was 
found in Ref. [21] for some positive-eigenvalue problems. 

IV. CONCLUDING REMARKS 



In this paper we have studied the behavior of the negative-energy eigenvalues of the Hamil- 
tonian (J5l) as functions of e. We find that the behavior of the negative-energy eigenvalues is 
completely different from that of the positive-energy eigenvalues. The positive-energy eigen- 
values remain real and positive on the infinite interval e > 0, and this is because the turning 
points never leave the Stokes wedges associated with this eigenvalue problem. However, the 
negative-energy eigenvalues occur in an infinite sequence of finite intervals, as in (f6|, and 
at the edges of these intervals the turning points leave the Stokes wedges. Moreover, while 
the positive-energy spectrum is entirely real, the negative-energy spectrum eventually be- 
comes complex when the energy becomes sufficiently negative except at the isolated values 
e = 2N-2. 

The smooth behavior of the positive energies and the choppy behavior of the negative 
energies of the Hamiltonian (|5|) bears a striking similarity to the behavior of the Gamma 
function T(z) for positive and negative z. The function T(z) is smooth and positive for all 
positive z. However, when z is negative, T(z) is only smooth on finite intervals of unit length, 
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FIG. 6: Comparison of the negative eigenvalues in the N = 3 region of e, 2 < e < 10, with the 
WKB asymptotic formula for these eigenvalues. The eigenvalues are plotted as diamonds (colored 
purple in the electronic version) and the WKB formula is plotted as solid curves (green in the 
electronic version) . At e = 4 the eigenvalues are exactly the negatives of the conventional positive 
eigenvalues of the sextic anharmonic-oscillator Hamiltonian H = p 2 + x 6 . When e/4 there are 
only a finite number of negative negative eigenvalues. 



and at the edges ot these intervals T(z) is singular. The finite intervals of e that we have 
found in this paper are also remarkably similar to the intervals of [1/K, 4K] (K = 1, 2, 3, . . .) 
found in Ref. [23], which presented a study of spontaneously broken classical VT symmetry, 
and the intervals found in Ref. [21] , which re-examined the work in Ref. [23] at the quantized 
level. 

To underscore the dramatic differences between the positive-energy and the negative- 
energy properties of the Hamiltonian (p)]), we examine the Hamiltonian at the classical level. 
In Figs. [7] and [8] we plot the classical trajectories in complex momentum space for e — 1, 3, 5. 
(We use momentum space here rather than coordinate space because in p space there are 
always two turning points while in x space the number of turning points varies with e [2~4"].) 
Note that the positive-energy classical trajectories are of uniform complexity and make 
simple loops around each of the turning points. In contrast, the negative-energy classical 
trajectories are extremely complicated, so complicated that to display them more clearly we 
include blow-ups of the trajectories in Fig. [Hj These classical trajectories strongly suggest 
that at the quantum level the structure of the negative-energy sector is more complicated 
and far richer than that of the positive-energy sector. 

We acknowledge conversations with J. Watkins. CMB is grateful to the Graduate School 
at the University of Heidelberg for its hospitality and he thanks the U.K. Leverhulme Foun- 
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FIG. 7: Comparison of the positive-energy and negative-energy classical momentum-space tra- 
jectories for the Hamiltonian ([5]) for e = 1, 3, 5. The positive-energy trajectories are relatively 
uncomplicated and make simple loops around the two turning points (red dots). The negative- 
energy trajectories are quite complicated (blow-ups of these trajectories are shown in Fig. pi). The 
contrast in complexity between the positive- and negative-energy trajectories is striking and sug- 
gests strongly that the negative-energy sector of the theory defined by the Hamiltonian (pi) is richer 
and more elaborate than the positive-energy sector. 
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FIG. 8: Blow-ups of the complicated portions of the negative-energy momentum space in Fig. [7j 



